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Abstract 

X-ray techniques have been used for more than a century to study 
the atomic and electronic structure in virtually any type of material. 
The advent of correlated electron systems, in particular complex ox¬ 
ides, brought about new scientific challenges and opportunities for the 
advancement of conventional x-ray methods. In this context, the need 
for new approaches capable of selectively sensing new forms of orders 
involving all degrees of freedom - charge, orbital, spin, and lattice - 
paved the way for the emergence and success of resonant x-ray scatter¬ 
ing, which has become an increasingly popular and powerful tool for 
the study of electronic ordering phenomena in solids. Here we review 
the recent resonant x-ray scattering breakthroughs in the copper ox¬ 
ide high-temperature superconductors, in particular regarding the phe¬ 
nomenon of charge order - a broken-symmetry state occurring when 
valence electrons self-organize into periodic structures. After a brief 
historical perspective on charge order, we outline the milestones in the 
development of resonant x-ray scattering, as well as the basic theo¬ 
retical formalism underlying its unique capabilities. The rest of the 
review will focus on the recent contributions of resonant scattering to 
the tremendous advancements in our description and understanding of 
charge order. To conclude, we propose a series of present and upcoming 
challenges, and discuss the future outlook for this technique. 
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1. INTRODUCTION AND BRIEF HISTORICAL SYNOPSIS 

Over the years, transition metal oxides have represented a traditional platform for strongly 
correlated electron physics, which has nowadays become a field of its own, encompassing 
several classes of compounds with one common denominator: the localized character of the 
low-energy electronic wavefunctions (with d or / orbital character) and the corresponding 
prominence of Coulomb interactions in driving the electronic properties of these materi¬ 
als. In conventional metals or semiconductors, the fermiology is essentially determined by 
the lowering of the total kinetic energy which becomes possible when a periodic potential 
supports the delocalization of the local orbitals into extended wavefunctions with a well- 
defined momentum and a correspondigly homogeneous distribution of the charge density. 
In correlated electron systems, the on-site Coulomb repulsion between two electrons in the 
same d or / orbital can overcome the kinetic energy part of the Hamiltonian, inducing the 
electronic system to find new ways to lower its total energy, often by spontaneous breaking 
of the native symmetries of the lattice (translational and/or point group symmetry). This 
tendency leads to the emergence of a rich variety of symmetry-broken electronic phases, and 
represents a distinctive trademark of strongly correlated systems, spanning across families 
of compounds otherwise very different from a chemical standpoint 0i- 

Within the extended class of correlated electron systems, copper oxides represent a 
unique breed due to the mixed character of the electronic bands, frustrated by the conflicting 
interplay between the 0-2p states - promoting itinerancy and electron hopping - and the Cu- 
3d states - which hinder charge fluctuations and are conducive to Mott-Hubbard physics and 
an insulating, antiferromagnetic (AF) ground state This delicate balance can be tuned 
and controlled by carrier doping, resulting in a phase diagram of astonishing richness and 
complexity, yet to be fully understood Antiferromagnetism has been long known as the 
ground state in charge-transfer insulating undoped copper oxides 0, while unconventional 
superconductivity was discovered in 1986 0. To date, several broken symmetries have been 
detected in the cuprates, which can be categorized into zero-momentum (Q — O) and finite- 
momentum {Q yf 0) orders, breaking rotational and translational symmetry, respectively. 
These different types of order can involve both the charge sector (nematic state for Q = 
0, and charge-density-wave for Q yf 0) and the spin sector (Q = 0 magnetic order, and 
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spin-density-wave or antiferromagnetism with Q^O), leading to several proposed forms of 
intertwined electronic orders |[^[6||14[). 


1.1. First observations of charge order in cuprates 

In this review we will focus on charge order, which is defined as an electronic phase break¬ 
ing translational symmetry via a self-organization of the electrons into periodic structures 
incompatible with the periodicity of the underlying lattice. The denomination of charge or¬ 
der, or equivalently charge-density-wave, has over time embraced a broad phenomenology. 
However, charge order was first discovered in the form of stripe order, whose real-space 
representation is depicted in Fig. In this ordered state, the doped holes (dark circles) 
are segregated into unidirectional structures which act as boundaries separating undoped 
domains characterized by AF order. Therefore, stripe order in its generalized meaning is an 
electronic ground state characterized by a combination of magnetic order and charge order 
with specific geometrical constraints on the ordering wavevectors. 

Historically, the first proof of stripe order came from neutron scattering, a momentum- 
resolved technique which had been extensively used in the early days of high-temperature 
superconductors (HTSCs) due to its excellent energy resolution and large magnetic cross 
section. The hrst neutron scattering studies of HTSCs focused on the doping and tempera¬ 
ture dependence of AF order in YBa2Cu306+a; (YBCO) (161 and La 2 _a:Sr 2 ;Cu 04 (LSCO) 
(17l, which was determined to be located at wavevectors Qaf = (±1/2, ±1/2, lQ Around 
the same time, the first experimental hint of stripe order was revealed by the independent 
neutron scattering measurements of Yoshizawa et al. (181 and Birgeneau et al. (191, who 
reported incommensurate AF modulations in underdoped Lai.92Sro.o8Cu04 (8% hole dop¬ 
ing) and Lai,89Sro.iiCu04 (11% hole doping), respectively, suggesting the presence of an 
electronic phase characterized by the quasi-static ordering of the doped holes. A subsequent 
study by Thurston et al. (201 revealed that the incommensurate AF order could only be 
found in the superconducting phase of LSCO [it was later shown to be present also in the 
insulating phase, albeit with different modulation vector (211]. The full momentum struc¬ 
ture and doping dependence of incommensurate spin order in LSCO was later uncovered 
by Cheong et al. ( |22[ ), and shown to manifest itself as a set of four magnetic peaks at 
Qaf-ic = (1/2 ± 5ic, 1/2, L) and Qaf-ic = (1/2,1/2 ± 5ic, L) (see diagram in Fig.[^), 
with 5ic representing the doping-dependent incommensurability. 

The experimental observations from neutron scattering spurred an intense theoretical 
activity, aimed at addressing two primary phenomenological findings: (i) the momentum 
structure of the incommensurate AF order and its relationship to the ordering of the doped 
holes; and (ii) the emergence of incommensurate AF order only within the superconducting 
phase, which pointed at an intimate interplay between short-range magnetism and super¬ 
conductivity in the cuprates. The short-range attractive interaction between segregated 
charges was initially proposed as a pairing mechanism for the superconducting state (231, 
while several numerical studies of the Hubbard model near half-filling were performed under 


AF: 
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High-temperature 
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YBCO; 

YBa2Cu306+a: 
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La2_a;Srj;Cu04 

LIMSCO: 

La2 —rc —y Hdy SrxCu044_6 


^Hereafter we will use reciprocal lattice units (r.l.u.) for wavevectors in momentum space. 
Reciprocal lattice units represent a notation where wavevectors are expressed as Q = {H,K,L), 
corresponding to Q = {H ■ 2-k (a, K ■ 217 fb, L ■ 2-k jc) in physical units (typically, A ^). Also, unless 
otherwise stated, we will refer the wavevectors to the undistorted unit cell, where the a and b axes 
(and correspondingly the reciprocal axes H and K) are parallel to the Cu-O bond directions with 
lattice parameters equal to the nearest Cu-Cu distance. 
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Figure 1 

Neutron scattering discovery of stripe order in 12% doped Lai 6 _ 2 ;Ndo. 4 Srj;Cu 204 . 

(a) Schematic representation of the stripe pattern: circles represent the Cu sites in the Cu02 
plane, with arrows denoting the Cu spins and gray circles indicating the location of the doped 
holes. (6) (HKO) projection of the reciprocal space showing the location of the measured charge 
and spin satellites (black circles) along the (010) and (110) directions, respectively; white circles 
represent Bragg reflections. {c,d) Elastic neutron scattering measurement of the charge order peak 
at (0,1.75,0) (c) and of the magnetic peaks at (0.5,0.375,0) and (0.5,0.625,0) (d), with scans 
along the arrows marked in (6). Readapted from Ref.|15| 


different conditions: using mean-field theory without (24 251 and with (261 charge fluctu¬ 
ations; within Hartree-Fock approximation (27 281; and by using exact diagonalization 
methods ( |29| ). All these studies pointed to a symmetry-broken ground state with short- 
ranged charge inhomogeneities organized into an ordered pattern of unidirectional charged 
arrays, segregated away from the magnetic domains and acting as boundaries for the latter 
(26). The stripe structure with holes condensed along rivers of charge-separating antifer¬ 
romagnetic domains was expected to manifest itself also as a periodic modulation of the 
charge, with a wavevector twice as long as in the case of the incommensurate AF peaks 
(24 261 and tied to the hole concentration (“a counting rule stating that the number of 
domain line unit cells is equal to the number of carriers.”, from Ref. 26). 
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Such a stripe state was found shortly thereafter by Tranquada et al. in the nickelate 
compound La 2 Ni 04+5 (301, as evidenced by the presence of charg^and magnetic satellite 
reflections in neutron scattering at the (collinear) wavevectors Qcharge = {H ± 2e, ± 2e, L) 

and Qspin = {H ± e, if ± e, L), respectively, reflecting the fundamental relationship between 
the spin and charge ordering incommensurability factors: ^charge = 2(5spin. These two charac¬ 
teristics (collinearity and two-fold proportionality of wavevectors) represent the basic con¬ 
ditions for stripe order, as opposed to other geometries (such as checkerboard). The search 
for static stripe order in the isostructural LSCO cuprate compound was complicated by the 
necessity for native fourfold symmetry in the Cu 02 planes to be broken in order to accommo¬ 
date a stripe state. A solution was eventually obtained with the addition of rare earth impu¬ 
rities (Nd and Eu, in place of La) distorting the lattice in the otherwise square-symmetric 
Cu02 planes and providing an inner strain field stabilizing stripe order. Consequently, 
in 1995 Tranquada et al. (151 used neutron scattering to detect charge and spin order 


peaks in the same sample of underdoped Nd-substituted LSCO, La 2 -a;-j,NdySra;Cu 204 +i 
(LNSCO) at Qcharge = (2, 2 ± 0.25,0) and Qspin = (0.5,0.5 ± 0.125, 0), as shown in Fig.[^ 
and d, respectively. This finding confirmed the wavevector relationship already proven in 
the nickelate compound, namely that 5charge = 0.25 = 2 • 0.125 = 25spin, as well as the equiv¬ 
alence between the incommensurability and the doping - with the former amounting to 
~ 1/8 = 0.125 at the ‘magical’ doping level of 12.5% - which is another signature of stripe 
order. This work represented the first direct evidence of stripe order in the cuprates and. 


together with closely following studies using again neutron (31 321, as well as x-ray (331 


scattering, initiated the whole experimental field of charge ordering phenomena in HTSCs. 

For more extended reviews of neutron scattering studies in the cuprates, we refer the 
reader to Refs. while for a comprehensive discussion of the theoretical works on charge 
and spin order we refer the reader to Refs.[^[7l|14[ 


1.2. Imaging charge order in real space 

Towards the end of the 90’s, a series of advancements in scanning tunnelling microscopy 
(STM) methods made it possible to obtain atomically-resolved maps with spectroscopic 
information on the energy-dependent local density of states (LDOS) via the measurement 
of the bias-dependent tip-sample differential tunnelling conductance as a function of the 
spatial coordinate r, dl/dV = g (r, V) oc LDOS {r,E = eV). The atomic resolution enabled 
by STM setups was key to the detection of density modulations with remarkably short 
correlation lengths, leading the way to a new era in the study of the nanoscale interplay 
between different electronic orders and of the resulting granular structure in HTSCs. 

The new spectroscopic imaging scanning tunnelling microscopy (SI-STM) capabilities 
were soon applied to the study of HTSCs, with the family of choice being that of Bi-based 
cuprates, owing to the presence of natural cleavage planes yielding extended, atomically-flat 
surfaces that turned out to be ideal for STM studies. The first compound to be studied was 
Bi2Sr2CaCu208+5 (Bi2212), a layered cuprate with a Bi0-Sr0-Cu02-Ca-Cu02-Sr0-Bi0 
stacking of crystallographic planes, characterized by a weak, van der Waals-type bonding 
between adjacent BiO layers. In 2002, 7 years after the original discovery of stripe order 
in the La-based compounds, the STM investigation of Bi2212 by Hoffman et al. (351 


^Neutrons can probe charge order indirectly, by detecting the associated distortion of the un¬ 
derlying lattice. 
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Figure 2 

Scanning tunnelling microscopy explorations of charge order in Bi-based cuprates. 

(a) Low-temperature topographic map of a cleaved surface of slightly overdoped (Tc =89 K) 
Bi2Sr2CaCu208+5 (Bi2212). (6) Differential tunnelling conductance (dl/dV) map (integrated 
between 1 and 12meV) between an applied magnetic field of 5 and OT, showing the appearance of 
period-4 density modulations within the magnetic vortex cores, and over the same spatial region 
of (a). Readapted from Ref. |35| (c) Momentum space representation of the periodic structures 
observed using STM on Bi2212, including the Fourier peaks from the atomic lattice (A) and the 
structural supermodulation (S), as well as the charge order peaks (Q) at Q ~ (±0.21,0) and 
Q ~ (0,±0.21), in reciprocal lattice units, (d) Energy-resolved Fourier-transformed conductance 
map along the (10) direction, showing no dependence of the charge order peak position (Q) on 
energy, thereby demonstrating the independence of this phenomenon from quasiparticle 
interference effects. Readapted from Ref. |36| (e-^-) Conductance maps in Na-CCOC at different 
bias voltages, and {h-j) corresponding 2D Fourier transforms, whose linecuts (k-m) along the 
Cu-0 bond directions show the charge order spatial frequencies and the lack of energy dispersion 
in the momentum-resolved structures. Readapted from Ref.|37| 


revealed the presence charge order in yet another cuprate family. The authors studied 
cleaved surfaces of slightly overdoped Bi2212 samples, whose topographic STM map is 
shown in Fig. HI- The corresponding differential conductance map over the same field of 
view, and upon application of a 5 T magnetic field, is depicted in Fig. if- In the field-induced 
vortex cores (darker spots), where the superconducting order is suppressed, four-fold, bi¬ 
directional modulations of the local density of states can be clearly imaged. This density 
modulation, with correlation lengths of the order of 30 A, provided experimental evidence in 
support of the expectation for a charge instability (381 in the doping range where the field- 
induced period-8 spin density modulation was discovered in LSCO shortly before ( |39[ ). The 
emergence of the periodic modulations within the vortex cores further suggested a direct 
competition between charge order and superconductivity, in analogy to the phenomenology 
in the La-based cuprates (40l. 


While similar evidence for charge order was found by Howald et al. (411 at low tem¬ 
perature and in the absence of a magnetic field, around the same time the discovery of 
quasi-particle interference (QPI) (42 I in the superconducting state and below an energy 
scale of the order of the superconducting gap underscored the importance of spectroscopic 


Riccardo Comin and Andrea Damascelli 




























measurements resolving the observed structure as a function of sample-to-tip bias voltage. 
The energy dependence of the charge order signal in Fourier space was later explored in 
the STM study of Vershinin et al. on Bi2212 (361 and Hanaguri et al. (371 on Na-doped 
Ca2Cu02Cl2 (Na-CCOC). Figure]^ shows the reciprocal space chart of Bi2212 - obtained 
by Fourier transforming the differential tunnelling conductance map - in the pseudogap 
state (T = 100 K), where QPI effects from superconducting quasiparticles are not active. 
The two-dimensional momentum structure of charge order can be also followed as a func¬ 
tion of energy (corresponding to the bias between the STM tip and the underlying surface), 
and Fig.[^ shows a series of linecuts across the charge ordering wavevector Q ~0.21 r.l.u., 
along the (HO) direction, as a function of energy. A non-dispersive charge order peak (Q) 
was found to be present across an extended energy range with almost constant amplitude, 
thereby proving its independence from the energy-dependent features due to quasi-particle 
interference. Charge order in highly underdoped (0.08 < p < 0.12) Na-CCOC, including 
non-superconducting samples, was discovered around the same period and shown to simi¬ 
larly exhibit an energy-independent momentum structure across a broad range of energies 
in the pseudogap state. Panels e-g and h-j in Fig.[^show, respectively, the two-dimensional 
real- and reciprocal- (Fourier-transformed) maps of the differential tunnelling conductance 
at 8, 24, and 48 mV bias. Distinctive features can be identified in the Fourier transform 
maps that correspond to near period-4 electronic modulations along the Cu-0 bond direc¬ 
tions. The energy-independence of these features appears evident from the linecuts shown 
in Fig.|^-m, demonstrating the presence of a broken symmetry driving static electronic 
modulations affecting all the electronic states over an extended energy range. 

The charge-ordered state in Bi2212 and Na-CCOC was later shown by Kohsaka et al. 

(431 to break down into bond-centered, locally unidirectional domains, as revealed by tun¬ 
nelling asymmetry imaging. These measurements also revealed the prominence of these 
electronic modulations at large energy scales, of the order of the pseudogap energy, a ten¬ 
dency confirmed in subsequent studies ( |44[|45[ ). Evidence for charge order in single-layered 
Bi-based compounds, Bi2yPbySr2-2La2Cu06+5 (Bi2201), was reported in 2008 by Wise et 
al. (461 over an extended range of dopings, demonstrating that the charge order wavevector 
evolution is compatible with a possible instability arising from the antinodal region of the 
Fermi surface. A similar doping dependence was also detected in Bi2212 (47 -491. In more 
recent years, a possible connection between charge order and the pseudogap state has been 
unveiled in Bi2212 using high-temperature STM to correlate the onset of the electronic 
signal from incipient stripes with the pseudogap temperature T* ( |50|||51[ ). These findings 
were followed by temperature-dependent measurements which show a competition between 
charge order and superconductivity near the Fermi energy (49l. Even more recently, it was 
shown that the zero-temperature density-wave order disappears in the vicinity of a putative 
quantum critical point located at a doping of p ~ 0.16 for Bi2201 (52 531 and p ~ 0.19 
for Bi2212 ( |54[ ). Throughout the years, STM and SI-STM provided a constant stream of 
scientific results and valuable insights on the existence, nature, symmetry, and nanoscale 
structure of charge order in underdoped cuprates; many of these studies have now been 
covered by several reviews (55 571. 


www.annualreviews.org • Resonant x-ray scattering studies of charge order in cuprates 
















RXS: Resonant 
X-ray scattering 


2. RESONANT X-RAY METHODS 
2.1. Resonant x-ray scattering in a nutshell 


X-rays have been long used to investigate the inner structure of matter, thanks to the in¬ 
teraction of light with the electronic clouds surrounding the atomic nuclei ( |58[ ). The role 
of photon energy has traditionally been secondary - with the exception of anomalous x-ray 
diffraction, a technique which relies on the strongly varying x-ray absorption near the ab¬ 
sorption edges of certain elements to simplify the phase retrieval problem in crystallography 
and provide a full reconstruction of the atomic positions in complex macromolecular sys¬ 
tems (1591). The first glimpses of resonant x-ray effects were revealed at the beginning of the 


1970’s, when de Bergevin and Brunei (601 demonstrated that x-rays are also sensitive to the 
distribution of electronic spins in magnetic materials by detecting antiferromagnetic Bragg 
reflections in NiO, thus confirming theoretical predictions by Platzman and Tzoar (61). 
Following these seminal findings, synchrotron-based x-ray magnetic scattering has been 
subsequently used as a powerful alternative to neutron scattering on several magnetically- 
ordered systems (62-67), while around the same time the foundations for a comprehensive 
theory of resonant x-ray scattering (RXS) were laid by Blume and coworkers (68-70). Since 
its inception during the 1980’s, resonant hard x-ray scattering has developed into an ex¬ 
tremely versatile tool (71 72) for the selective study of ordering phenomena involving the 
charge, spin, orbital, and lattice degrees of freedom, often offering a unique perspective on 
their respective interplay. 

Historically, hard x-ray methods have anticipated soft x-ray methods due to the (huj)^ 
dependence of the x-ray attenuation length on photon energy {to is the angular frequency 
of the radiation field), which eliminates the need for vacuum-based experimental chambers 
for energies hu > 5 ke\jf] Soft x-ray absorption (XAS) was first pioneered using conventional 
x-ray sources and in a transmission geometry on thin films of rare earth metals and oxides, 
known for the large absorption cross-section and rich multiplet structure characterizing the 
M 4,5 edges (3d —4/ transitions) ( |75[|76[ ). Important advancements in XAS came with the 
use of synchrotron facilities ( |77[ ), which produced a higher x-ray flux and further enabled 
the control of linear and circular x-ray polarization, making it possible to detect magnetic 
x-ray dichroism effects (78). In the same years, improved detection schemes were being 
progressively adopted, such as partial electron yield (791, total electron yield ( |80[|81[ ), and 
fluorescence yield (|82l|83l). 


The successes of resonant hard x-ray scattering spurred the development of soft x-ray 
scattering instruments and methods, which presented clear benefits - access to the electronic 
states and corresponding degrees of freedom controlling the electronic structure in transi¬ 
tion metal and rare earth oxides, greater sensitivity to surface and interface effects, native 
control of incoming light polarization - but also inherent complications - limited accessible 
reciprocal space (and a consequent lower boundary of ~ 3 Aon the smallest measurable peri- 


®For example, the attenuation length for x-rays propagating in atmosphere at lOkeV (IkeV) is 
~3m (3mm) and, as a result, the loss in x-ray flux over 1 m (representing the typical dimension 
of a scattering diffractometer) is ~ 30% (100%). Due to the continuous, power-law dependence of 
x-ray absorption on the photon energy, no precise cutoff can be defined for soft x-rays. However, 
conventional Cr-Ka x-ray sources exist that operate in air at ~ 5.4 keV, albeit requiring reduced 
source-to-sample and sample-to-detector distances. In addition, we note that hard x-rays beamlines 
still require vacuum conditions along the long pipes transporting the photons from the storage ring 
to the experimental chamber and from the chamber to the photon detectors. 
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Figure 3 

The first resonant soft x-ray scattering studies of charge order in cuprates. (a,b) 
Reciprocal space mapping of the (HOL) plane in the non-superconducting cuprate compound 
Sri4Cu2404i measured using soft x-ray scattering (a) off-resonance (526eV) and (b) at the main 
absorption peak for doped 0-2p holes (528.6eV, corresponding to the Is —> 2p transition). At 
resonance (b), a superlattice peak is revealed at Q = (0, 0, 0.2), signaling the crystallization of 0-2p 
holes into a static and periodic pattern (Wigner crystal). Readapted from Ref. |73[ (c) Resonant 
soft x-ray scattering measurements of the stripe order peak of Lai, 875 Bao.i 25 Cu 04 [located at 
Q~ (0, 25, 0,1.5)1 across the O-K edge (Is —> 2p, left), with highlighted the x-ray transitions into 
the mobile carrier states (mobile carrier peak, MCP) and upper Hubbard band (UHB). (d) Same 
as in (c), but at the Cu-Ls (2p —> 3d, right) absorption edge. The scattering peak intensity (red 
line and markers) is resonantly enhanced in the vicinity of the features of the absorption spectra 
(green line) corresponding to electronic transitions into the O and Cu sites in the Cu02 planes. 
Readapted from Ref.|74| 


odicity, at ~2keV), the incompatibility with crystal-based energy or polarization analyzers, 
the need for high vacuum environment. The rise of resonant soft x-ray scattering required 
special diffractometers capable of simultaneous, in-vacuum control of the sample and de¬ 
tector angles. The first concept of a two-circle diffractometer (with independent sample 
and detector rotations in the scattering plane) was pioneered already in the late 1980’s ( |84[ ) 
and further advancements were realized throughout the 1990’s and early 2000’s at different 
facilities including Brookhaven National Lab (Ref. |85| and later Ref. 861, LURE (Ref. |87[ ), 
the Synchrotron Radiation Center (Ref]88[), the European Synchrotron Radiation Facility 


891, the Daresbury Laboratory (Refs. |9d||92[) , BESSY (Ref. |93[ ), the Swiss Light Source 
941, and the Canadian Light Source (Ref.|95[) 


(Ref. 

(Ref. 941, and the Canadian Light Source (Ref. |95[ ). Motivated by the new insights produced 
by early resonant inelastic scattering studies at the Cu-K edge (~8.9keV) (|96||99|, the first 


resonant scattering measurements on cuprates in the soft x-ray regime were performed in 
2002 by Abbamonte et al. ( |100[ ) on thin hlms of oxygen-doped La2Cu04, at the O-K edge 
(Is —>■ 2p, fia;~530 —550eV) and Cu-L3,2 edge (2p —>■ 3d, fia;~920 —960eV). The sensitivity 
of these two absorption edges to the doped holes in cuprates had been demonstrated previ¬ 
ously ( |8T1 ioihio 3P , in particular through the identification of a mobile carrier peak (MCP) 
structure at fiu;~528eV, i.e. below the onset of the main O-K edge, and representing tran¬ 
sitions onto the doped hole electronic states with large spectral weight on the 0-2p orbitals 
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within the CUO 2 planes. The study by Abbamonte et al. explored momentum space in the 
ranges (0, 0,0.21) — (0,0,1.21) (orthogonal to the CUO 2 planes) and (0,0, 0.6) — (0.3,0,0.6) 
(parallel to the CUO 2 planes) and the evolution of the momentum-dependent interference 
fringes across the resonances was determined to be indicative of a rounding of the carrier 
density near the film-substrate interface. Most remarkably, the authors revealed the extent 
of resonant enhancement at these absorption edges in the cuprates, which amounts to a 
single doped hole in the CUO 2 planes scattering as strongly as 82 electronic charges, with 
a resulting magnification of the experimental signal equal to the squared of the scattering 
amplitude, or 82^ > 10^. This can be regarded as one of the key figures of merit for resonant 
x-ray methods, and represents the foundational mechanism underlying the success of RXS 
in detecting weak ordering phenomena (such as charge order) in the cuprates. 

Later, Abbamonte et al. ( |73| ) used RXS to discover and characterize the ordering of 
doped holes in the non-superconducting cuprate compound Sri4Cu2404i. A peak in re¬ 
ciprocal space, representing the signature of a well-defined, periodic modulation of the 
electronic density along the reciprocal L axis and with wavevector Q = (0,0,0.2), could 
only be measured at the O-K pre-peak resonance at 528.6 eV (Fig.j^) while no signal was 
detected once the photon energy was tuned off-resonance by only a few eV (Fig.[^). This 
work represents the very first RXS evidence of charge order in a copper-oxide compound, 
corresponding to the crystallization into a Wigner crystal phase, manifested as a modu¬ 
lation of the electronic charge triggered by electron-electron Coulomb interactions. This 
observation demonstrates the power of resonant scattering methods as opposed to conven¬ 
tional diffraction techniques, which, being non-resonant, are largely insensitive to subtle 
electronic ordering phenomena not involving the lattice degrees of freedom. 

The hrst RXS study of charge order in one of the superconducting cuprate families came 
along in 2005, when Abbamonte et al. (741 investigated stripe order in LBCO. While the 
momentum structure of the stripe-ordered phases in this compound had been previously 
studied using neutron scattering (321, no direct information was yet available on the role 
and involvement of the electronic degrees of freedom, due to the fact that neutron scattering 
predominantly measures periodic distortions in the lattice. Abbamonte et al. found a peak 
in reciprocal space at (0.25,0,L), consistent with previous observations (321 and with a 
weak dependence on the out-of-plane momentum L, thus conhrming the prominent two- 
dimensional nature of the charge ordered state. Once again, this study reaffirmed the 
fundamental role played by the resonant enhancement: Figure]^ and d show the intensity 
of the charge order peak (red markers) overlaid onto the absorption profiles (green lines) for 
the O-K and Cu-La edges, respectively. The intensity for charge scattering is nonzero, i.e. 
detectable on top of the fluorescent background and above noise level, only near the Oxygen 
prepeak (MCP) and the Cu excitonic resonance. This work had profound implications 
for the understanding of charge order in cuprates, since it clarified that stripe order is 
predominantly configured as a modulation of the electronic density followed by a distortion 
of the lattice, which however represents a secondary effect, as demonstrated by the absence 
of any diffracted intensity away from the electronic resonances (in the soft x-ray range). 

In the following sections we will provide a brief description of the theory of resonant 
scattering and of the experimental scheme. For a more comprehensive treatise on these 
topics, we refer the reader to Refs. [68l[7Tl[7^ pMl[T05| and references therein. 
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2.2. Theory of resonant scattering 


Resonant x-ray scattering is a photon in - photon out technique, where photons get scattered 
from a material due to the interaction with the electronic clouds. For radiation-matter 
scattering to occur, the interaction Hamiltonian has to contain operator combinations of 
the kind a^, (q) al (q — Q), where al (q) [oi, (q)] is the operator creating [annihilating] a 
photon with wavevector q, polarization state v, and frequency uj = cjqj. The effective 
nonrelativistic interaction Hamiltonian can be derived from the full electron-matter minimal 
coupling Hamiltonian, and reads: 


Htot = 


E, {i [P^ - f A(r„ t)]^ + H(r„ t)} + E, 


j^k |rj-rfc|2 


+ 


— -^^el + -^^EM H- 'y ^) ' Pj “1“ 

'-V-' rrieC^^j 


2meC^ 


• ^ 1 


( 1 ) 


Ho 


where e and m are the fundamental electronic charge and mass, pj and represent 
the momentum and position coordinates of the j-th electron respectively, and V (r, t) and 
e^/jr — r'j are the lattice potential and the Coulomb interaction terms, respectively. A(r, t) 
represents the vector potential, Jfei = Ej f) + Ej^fc |r -rfcF Hamilto¬ 
nian of the electronic system, while i?EM = Eq [“t (^l) (q) + 1/2] is the Hamiltonian 

of the electromagnetic (EM) field alone. 

The interaction operators and which are respectively linear and quadratic in 

the vector potential, couple the electromagnetic field and the electronic degrees of freedom. 
At this point, we can use as basis set for the light-matter quantum system the states 
|^m)= X where represents the electronic part of the wavefunction 

(with eigenvalues em, and m labeling a generic set of quantum numbers) while 
indicates a photon state with photon occupation hq,v = ...}, corresponding 

to having photons with wavevector and polarization (qi,izi), nq 2 ,v 2 photons with 

(qzjJ^z), and so on. In this notation, M labels the global set of quantum numbers, i.e. 
M = {m, q, u}. Note that due to the radiation-matter interaction, the states |^'m) are not 
eigenstates of the system, but they can be used as basis set in a perturbative scheme, in 
which case we define the unperturbed (i.e., with the interaction terms turned off) energy 
spectrnm as Em =Cm + Eq v {hci,+ 1/2). Within this framework, a scattering process 
is defined as a transition from an initial photon state \4'i)EM ~ l■■■)l^qin,‘^in)l^qout,!^out)l■ ■) 
to a final photon state \4>f)EM ~ l■■■)l^qin,‘'in “ l)l’^qout,i'out + 1)|...)) corresponding to the 
annihilation of an incoming photon with wavevector and polarization (qin,J^in) and the 
concomitant creation of an outgoing photon (qout, t'out). 

Central in the theory of elastic scattering is the calculation of the probability of tran¬ 
sition from a state |5'i)= |'!/’Gs)ei x l</i)EM 1° state |^'/)= |'!/'Gs)ei x |<//)bmi where the 
photon states |<(>i)EM I/’/)em given in the previous paragraph, and where we 

further assume that the electronic part of the initial and final state is in the ground state 
|'*/’Gs)ei- The transition probability Wi^f between the initial and final quantum states can 
be calculated using the generalized Fermi’s golden rule ( |106[ ): 

Wi^f=27T\{tS>^\T\tS>f)f5(Ef-E.) (2) 


where the T-matrix is defined as follows: 

1 


T = Hint + Hu 


' Ei — Ho + irj 


Hint + Hu 


-Hi 


Ei-Ho + iv E,-Ho+irj 


Hint + ... (3) 
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Figure 4 

Resonant processes and scattering geometry in RXS. (a) In non-resonant scattering the 
excitation process does not involve intermediate states, while resonant scattering occurs whenever 
the incident photon energy is tuned to promote an electronic transition from the ground state '^gs 
to an intermediate state The subsequent radiative recombination of the excited electron with 
the core hole results in the creation of an outgoing (scattered) photon, (ft) The different photon 
energy dependence of resonant and non-resonant processes, showing the enhancement occurring in 
the resonant channel near an electronic transition with energy AE. (c) Schematics of a 
conventional diffractometer, with kinematics for the scattering/diffraction process outlined in (d). 


Here the first operator on the right-hand side represents the first-order perturbation term, 
the second operator represents the second-order perturbation term, and so on. In Eq.[^ 
Hint is the interaction operator and Hq is the unperturbed Hamiltonian - in our case 
Hint =^^int + ^0 = Hgi + H^M, respectively. In scattering, we require operator 

combinations of the kind a^a inside the T-matrix, which originate from interaction terms 
that are quadratic in the vector potential A, since the latter can be expressed in second- 
quantized notation as A(r, t) oc ■ [exp(iq • r — iuit) al (q) + h.c.] (e^ denoting the 

polarization vector of the polarization state u). These combinations are generated by using 
the quadratic interaction operator in the first-order term of Eq.|^ and by the linear 

interaction operator H]^ in the second-order term of Eq. 3 The corresponding first 
(‘ 2 '\ 

and second order perturbative transition probabilities can then be obtained from 
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= 271 


wf2f = 271 
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rrieC 


E, 


Ei — Em + *r M 


(5) 


where for convenience we have dropped the 5-function (enforcing conservation of the total 
energy in the scattering process) originally present in Eq.[^ 

In Eq.[^ represents a generic (excited) quantum state of the light-matter system, 
with corresponding energy Em and lifetime Ti/Tm- The squared vector potential can then 
be expanded as: 

A^(r, t) oc |^gdq >—+ /i.c.j x aj,, (q') -|- /i.c.j 

« (e-'in • (qout) • (qin)] 

oc (e„j„ • Ei.„„t) X {q_out) (qin) (6) 


where and Sj/out represent the polarization vector of the incoming and scattered photons, 
and in the last step of Eq.j^we assumed the scattering process to be elastic, i.e. tuin =iUout- 
At this point, by using the previous definitions for the initial and final states |5'i) 
and Ivlt/), together with the fact that Ei = ecs + [’^qin^^qin + 1/2] and Em = tm + 
[(^qin - 1) ^^qin + 1/2], and considering that (qout) aj,.^ (q^) |^i)EM ^ 
rewrite Eqs.[^and[^as follows: 


w 


( 1 ) 


oc 


w 


( 2 ) 


OC 


(V^Gsiy^ e |' 0 Gs)| oc KV^GsI p (Q) |' 0 Gs)|^ 

tGS - tm huj irm 

m j,k 


2 


(7) 

( 8 ) 


where Q = qin — qout is the momentum transferred by the photon field to the sample, and 
p (Q) is the Fourier transform of the electron density operator p (r) = Ej'^ (*■ “ *’i)- 

In the x-ray regime, the transition channel represented by Eq.|^ involves the excitation 
of a high-energy many-body state with a core hole (tpm)- In a more intuitive perspective, 
this is equivalent to the excitation of a core electron into an intermediate state through 
absorption of the first photon, followed by re-emission of a (scattered) photon once the core 
hole is filled back. On the other hand, in the first-order perturbative term expressed by 
Eq.0 the scattering process is instantaneous as it does not involve the excitation of an 
intermediate state. As a consequence, the first-order mechanism (known as Thomson scat¬ 
tering), proportional to the squared amplitude of the total electronic density in the ground 
state, is non-resonant and controls the signal in conventional x-ray diffraction (XRD). The 
second-order process is instead resonant, and is therefore associated to RXS. 

Equation]^ is a particular version of the Kramers-Heisenberg formula, which represents 
the general solution to the problem of a photon scattering from an electron. It can be 
further simplified under the assumption that the x-ray excitation process is local, which 
implies that: (i) the local orbitals (at the lattice site n) can be used for the electronic basis 
set: ■!/>m(r) —>• = Xi(r “ R-n); (ii) that all matrix elements oc Sm,n, 


XRD: X-ray 
Diffraction 
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i.e. that they vanish for orbitals belonging to different lattice sites; and (iii) that, due to 
the localization of the initial (core) electron around the position R„ of its parent atom, the 
phase due to the photon held can be approximated as At this stage, it is 

convenient to introduce a new quantity, the form factor fpq, which is a photon energy- and 
site-dependent complex tensor dehned as follows: 


fpq^ (fi^) 






hoj - (4"' - -I- iFii 


(9) 


Here Xi"'' ^nd xi"^ represent the initial and intermediate single-particle electronic states at 
site R„ (with energies and respectively) involved in the light-induced transition 
i ^ 1. Fii is the inverse lifetime {h/ru) of the intermediate state with an electron in 
and a hole in • 

The resonant scattering cross-section, through the form factor fpq^ (fiw), can be shown to 
bear a close connection to the x-ray absorption (XAS), which is a hrst order process in the 
radiation-matter interaction Hamiltonian: 


(huj) oc - ^ X Im ■ fpp^ (^^)] (10) 

7^^® (Q,fitu) oc |^^^(e,.i„)p • l^^fpq^ (huj) e '^'• (e^o„t)g| 

= |5Zp,j(^'^i")p ■ (^‘^) ■ • (11) 


In Eq. 11 Fpq represents the scattering tensor, which is not a local quantity (does not 


depend on the lattice position R„) and is more directly related to the physical observable 
in RXS experiments (7^^®). Moreover, from the above equations it follows that XAS only 
depends on the incoming light polarization whereas the RXS signal depends on the 
outgoing light polarization as well. 

The difference between resonant and non-resonant scattering is schematized in Fig.|^a). 
The mechanism corresponding to XRD involves a single step, in virtue of its hrst-order 
nature; conversely RXS, being a second-order transition, proceeds in two stages involving an 
intermediate state. This clearly reflects in the very different photon energy (hv) dependence 
of the two channels [see Fig.|^b)]: whereas XRD is nearly energy-independent (red dashed 
curve), the cross-section for RXS is strongly peaked around the energy of the electronic 
transition (blue curve), where the experimental signal undergoes a strong enhancement 
(while decaying to zero away from the resonance). This usually occurs in correspondence of 
an absorption edge, i.e. when electronic transition from a deeply bound core state into the 
valence band (and beyond into the continuum) take place. As a consequence, RXS gains 
a strong sensitivity [as large as 82-fold in the cuprates 174 1001] to partial modulations of 
the charge density involving a single electronic band, whereas the XRD signal reflects the 
total electronic density (see Eq.[^ and therefore suffers from a weak sensitivity to spatial 
variations of the latter, unless they are accompanied by a distortion of the lattice, which 
would involve all the electrons (core and valence). 


2.3. The experimental scheme 

In an actual scattering or diffraction experiment, a monochromatic x-ray beam with 
wavevector qin, photon energy hcjin = c ■ qin and polarization impinges on a sample 
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and a scattered photon will be detected along the direction of the wavevector qout using an 
energy-integrating photon detector (PD) or an energy-resolving spectrometer (see Fig. &)• 
At the end of the process, a net momentum and energy have been transferred to the sample, 
which can be derived from the corresponding conservation laws: 

hi/in = huout + AE (12) 

qin = qout + Q. (13) 

In the case of elastic scattering hvin = huout and there is no energy transfer with the sample 
(Ai? = 0), whereas the case AEj^O defines inelastic scattering events. Strictly speaking, 
elastic scattering probes the static component of the charge and magnetization density 
occurring in the system under study, whereas inelastic scattering is sensitive to dynamical 
processes and low-energy excitations. However, due to the spectrometer-characteristic finite 
energy resolution SE (which in the soft x-ray regime ranges between 30 meV and 1 eV, 
while hard x-ray spectrometers reach down to ImeV), purely elastic scattering cannot be 
experimentally accessed, and it is more appropriate to use the term quasi-elastic scattering, 
which probes a regime which is static up to a timescale r ~ h/5E. From a more practical 
standpoint, the energy-integrated measurement in most cases yields a reliable representation 
of the momentum structure of the ordered state, due to fact that the inelastic part of 
the spectra usually evolves very smoothly and can be discarded as background in RXS, 
especially if it exhibits a different temperature dependence with respect to the zero energy- 
loss feature (see also discussion of Fig.j^. 

Hereafter, we will focus on the momentum structure of RXS measurements and, un¬ 
less otherwise specified, will assume the use of energy-integrated mode. From Eq. |13| 
and using hi/in = hvout, the magnitude of the exchanged momentum can be expressed 
as Q = 2qin X sin(6^sc/2) [0sc is the scattering angle, see Fig.j^c)]. Projecting Q into the 
plane defining the sample surface then yields the in-plane (Qy) and out-of-plane (Q±) com¬ 
ponents of the transferred momentum, which will be often referenced in the rest of this work 
[see Fig.j^d)]. The realization of this geometry is illustrated in Fig.|^c), and described in 
greater detail in Ref. |95| In general the photon detector moves on a single-circle, that is, a 
single angular goniometer (the corresponding variable is often denominated 26 and corre¬ 
sponds to the scattering angle 0sc). The sample stage usually involves translational motion 
(xyz) and various rotations, whose number defines the type of diffractometer. Two-circle 
diffractometers represent the most common choice for soft x-ray experiments, featuring a 
single angular motion (with angular variable denominated 6) for the sample (F'^ circle), 
whose axis is perpendicular to the scattering plane (the one spanned by the vectors qin and 
qout), in addition to the detector rotation (26), which represents the 2"*^ circle. The first 
generation of soft x-ray diffractometers also includes two additional rotational degrees of 
freedom (denominated x and (j), see again Fig.|^), allowing a fine alignment of the sample 
axes with respect to the scattering plane (however, x and (j> typically cover a limited angular 
range of —5° to 5°). New designs are being developed using different geometries to extend 
the angular range and control for the sample orientation. Diffractometers with more circles 
(up to 6 - 3 for the sample, 3 for the detector) covering an ample angular range are routinely 
used at higher photon energies (hard x-rays, hw > 7keV), where the added complication 
of the all-vacuum environment required for soft x-rays is lifted. Consequently, hard x-ray 
diffractometers can typically access a wider portion of reciprocal space. 

Typically, the experimental signal is comprised of both resonant and non-resonant con¬ 
tributions, and in the two possible regimes or ^ one ends 
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up probing different phenomena (see again Fig. &)• In the first case, where non-resonant 
processes are dominant, ail electrons contribute equally to the measured signal, which is 
therefore simply proportional to the atomic number Z. The diffraction signal will be dom¬ 
inated by the core electrons, which usually outnumber the valence ones (ricore Uvaience), 
with the exception of lighter elements which as a result are not probed very effectively in 
XRD. In addition, since core states are very tightly bound to their parent nucleus, in a 
conventional diffraction experiment one mainly probes the ionic lattice in reciprocal space, 
which is why XRD is used primarily for structural studies. In the second case, the scattering 
process has a strong enhancement in correspondence of a very specific electronic transition. 
As a result the signal bears the signature of the electronic distribution (in reciprocal space) 
of the final state of such transition. This characteristic of resonant scattering allows it to be 
not only element-specific (whenever the absorption edges of different chemical species are 
spaced sufficiently apart in photon energy), but also orbital-selective. This unique capabil¬ 
ity of RXS has been established and employed in many different systems. Charge-ordering 
in cuprates (74) and cobaltates (1071, and orbital-ordering in the manganites ( ^|108|[l09 1, 
are among the most spectacular case studies. 


3. CHARGE ORDER IN CUPRATES - A NEVER-ENDING JOURNEY 
3.1. A resurging phenomenology: charge-density-waves in YBCO 

Fifteen years after the original discovery of stripe order by Tranquada et al., charge order 
had been observed in the doping region around 12% hole doping directly in the 3 families 


of La-, Bi- and Ca-based cuprates: LNSCO ( 15||3ip3l40yil0 1, LBCO |32 

74l|110||lll|, and 

LESCO |ll2l 113| using neutron and x-ray scattering (in the case of Refs 

^^|ll3 resonant 

x-ray scattering was used, in particular); Bi2212 (35 -371141 43||50||51 114 

■nolTNa-ccoc 


(37 431, and Bi2201 (461 using STM. However, no clear evidence of charge order had been 


LESCO: 


found in the YBCO compounds, where the possibility of introducing doped carriers via the 
La 2 ” Eu Sr CU 2 O factional filling of the Cu-O chain layer helps reducing the electronic inhomogeneity in the 

The first indication of electronic order in YBCO came in 2007 thanks to a series of 
pioneering measurements at high magnetic fields (up to 62 T) and low temperatures (down 
to 1.5 K) revealing the presence of quantum oscillations in the Hall resistance in the 


normal state of underdoped YBa2Cu306.5i (p = 0.10, Tc = 57K) samples (122 1231. The 
frequency of these oscillations vs. the inverse magnetic field provided evidence for the 
emergence of small Fermi pockets at high fields ( |122| , arguably of electron-like nature in 
light of the negative sign for the Hall coefficient Rh ( |123| . These findings indicated a change 
in the Fermi surface topology from the large hole-like Fermi ‘barrels’ in the overdoped regime 
( 124|125 1 to smaller pockets in the underdoped regime, thereby suggesting “a reconstruction 
of the Fermi surface caused by the onset of a density-wave phase, as is thought to occur in the 
electron-doped copper oxides near the onset of antiferromagnetic order” (|123|). The striking 


similarity between the Hall and Seebeck coefficients of YBCO and those of LESCO (126 


1281 made a strong case for a form of order in YBCO akin to the stripe order in LESCO. 


The first direct evidence that the Fermi surface reconstruction is due to charge order was 
provided by high-held nuclear magnetic resonance (NMR) measurements on YBa2Cu306.54 
(p = 0.104, Tc = 57K), where a splitting in the ®®Cu(2) lines was observed at 28.5T and 
below 50 K, signaling a change in the quadrupole frequency which was ascribed to a periodic 


variation in the charge density at planar Cu sites or at the oxygen sites bridging them (1291. 
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Figure 5 

Resonant x-ray scattering discovery of charge-density-waves in (Nd,Y)Ba2Cu306+a;* 

(a) Resonant inelastic x-ray scattering (RIXS) measurements of underdoped Ndi. 2 Bai. 8 Cu 307 
(T’c = 65K), as a function of energy loss and momentum [along (ffOO)], showing the emergence of a 
quasi-elastic peak around ~—0.31 r.I.u.; inset: three-dimensional view of the lattice structure of 

YBa2Cu306+cc- Momentum dependence of the quasi-elastic (6) and of the energy-integrated 

RXS intensity (c) for a series of temperatures across the superconducting transition Tc. {d,e) 
Temperature evolution of the charge order peak intensity (d) and full-width-at-half-maximum (e), 
showing a cusp at Tc, thereby providing evidence of competition between charge order and 
superconductivity. Readapted from Ref. El] 


The simplest scenario seemingly compatible with the NMR results was a unidirectional 
density modulation with 4a (four unit cells) periodicity. Furthermore, the observation of 
the charge-order-induced NMR splitting for strong fields perpendicular but not parallel to 
the conducting Cu 02 planes provided evidence of a competition between superconductivity 
and charge order. 

The quantum oscillations and NMR results seemed to point to the necessity of com¬ 
pletely suppressing superconductivity to observe the emergence of a seemingly competing 
charge ordered state. Following early RXS explorations by Hawthorn et al. (1301, the first 
direct observation of charge density wave in reciprocal space, and at zero magnetic field, was 
obtained in 2012 by Ghiringhelli et al. (121 \ using energy-resolved and energy-integrated 
RXS on a series of YBCO doping levels around 12%. The central experimental data uncov¬ 
ering the momentum location, and therefore the periodicity, of charge order in YBCO are 
shown in Fig. §i, and consist of a series of scans of the x-ray energy loss (horizontal axis) 
for different values of the planar projection (H) of the momentum Q = (JT, iT, L) along the 
(HOL) direction in reciprocal space. The energy structure of the resonant inelastic x-ray 
scattering spectra bears three main contributions: (i) a quasi-elastic line at zero energy 
loss; (ii) a low-energy peak/shoulder representing spin excitations; and (iii) intra-band (dd) 
particle-hole excitations. While the spectral weight of the spin and dd excitations is nearly 
momentum-independent, a clear enhancement of the quasi-elastic line can be seen around a 
planar momentum 0.31 r.l.u., which reveals the presence of a periodic modulation of 

the electronic density - representing the momentum structure of a zero-field ordered state 
possibly connected to the one previously identified with NMR [recent high-held diffraction 
measurements, which will be discussed later, helped clarify the nature of this connection 
(131l]. The RXS data further revealed charge modulations to be present along both pla¬ 
nar crystallographic axes, albeit with different amplitudes (as discussed more extensively 
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in Ref. 1321. The same momentum-space peak could be equally well identified following 


the quasi-elastic component in energy resolved spectra (Fig.[^) and the energy-integrated 
RXS (Fig.j^), which confirms: (i) the importance of the resonant enhancement at the Cu-I /3 
edge; (ii) only the quasi-elastic scattering contributes to a pronounced momentum-resolved 
structure, as previously discussed. 

The detailed temperature dependence of energy-resolved and energy-integrated RXS 
momentum scans in Figs. |5> and c points to an onset temperature T ~ 150 K and, most 
importantly, to a partial suppression of the charge order peak below the superconducting 
transition temperature Tc. This finding provides direct evidence of a competition between 
superconductivity and charge order, and is further substantiated by the temperature evolu¬ 
tion of the charge order peak intensity and full-width-at-half-maximum (FWHM, inversely 
proportional to the correlation lengths) reported in Fig. ii and e for the case of energy- 
resolved and energy-integrated RXS, respectively. The intensity and FWHM data indicate 
a clear weakening of the charge ordered state both in its amplitude and spatial correlations 
at the emergence of the superconducting state. Around the same time, the wavevector of 
the charge-density-wave in YBCO, and its competition with the superconducting state were 
measured by non-resonant hard x-ray diffraction experiments, which revealed an enhance¬ 
ment of the charge order amplitude below Tc when superconductivity was actively weakened 


by applying magnetic fields up to 17 T (133l. Field-induced enhancement of charge order 


was similarly reported in LBCO, using hard x-ray diffraction (134l. 


The discovery and identification of the charge order state in YBCO represented a break¬ 
through in the field of copper-oxide superconductors, as it suggested that the charge ordered 
state could be a defining instability of the CUO 2 planes, and strongly contributed to re¬ 
vitalizing this research direction. Stimulated by the possibility to provide a unifying view 
of charge order in the cuprates, intense efforts at the experimental and theoretical level 
were put forth in order to identify a common thread across the different manifestations of 
charge order that emerged over the years. A series of studies analyzed the photon energy 
dependence of the charge order signal from LESCO ( |112[ ), LNSCO ( |135[ ), and YBCO ( |136[ ), 
which is inherited from the photon energy dependence of the complex form factor / (fitj) 
(see Eqs.j^and 111. The imaginary part of the latter is directly related to the x-ray ab¬ 
sorption signal, and the real part can be extracted by a Kramers-Kronig transformation. 
Eor a transition into a single C\x-‘id^ 2 _y 2 hole (which is the case at the Cu-Z /3 edge in 
the cuprates) there is no multiplet structure, and the (site-dependent) form factor can be 
approximated by a simple Lorentzian lineshape as ~ An{hio — e„ + iF)”^, so that the 
RXS intensity can be expressed as: 


(Q, hut) oc {huj) 

‘ ^ n. 


E 


An 


iQ* Rn 


n. (hoj — £n + iF) 


(14) 


where the polarization dependence has been neglected for simplicity. Eq. |14| shows that there 
are three variables controlling the RXS signal which can vary spatially: (i) the transition am¬ 
plitude An', (ii) the lattice position R,,; and (iii) the transition energy Correspondingly, 
Achkar et al. evaluated the impact on the RXS signal of assuming a periodic modulation 
for (see Fig.[^): (i) the valence modulation; (ii) the lattice displacements; (iii) the transi¬ 
tion energy shifts (this term was introduced in Refs 135||136 1 . The comparison between the 
experimental RXS peak amplitude and the model calculations vs. photon energy is shown 
in Fig.|^ for LNSCO and in Fig.[^ for YBCO. In both cases, it was determined that the 
RXS signal is predominantly controlled by the spatial variation of the x-ray transition en- 
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Figure 6 

Photon energy-dependent resonant scattering from charge order in cuprates, (a) 

Schematic representation of the possible microscopic contributions to the RXS signal: valence 
modulation (top); lattice displacement (middle); transition energy shift (bottom). (6) Comparison 
of the photon energy-dependent RXS intensity from the stripe order peak in LBCO to calculations 
from the models introduced in (a), showing that the resonant scattering signal is predominantly 
controlled by spatial variations in the x-ray transition energies. Readapted from Ref. |135l (c) 
Photon energy-dependent RXS signal from charge order in YBCO for horizontal (tt) and vertical 
(c) incoming polarization (top) and comparison between experiment and the prediction of the 
energy-shift model for polarization along the a axis. Readapted from Ref.|136| 


ergy shifts, while in the earlier work on LESCO the RXS lineshape was explained with the 
nonlinear increase of the charge-carrier peak as a function of doping p, due to a reduction 
of correlation effects (U) at higher doping concentrations. More recently, an alternative 
theoretical framework has been devised that explains the photon energy-dependent RXS 
signal by accounting for the delocalized character of intermediate states ( |137[ ). 

A closely following series of studies provided a complete mapping of the detailed char¬ 
acteristics of this phenomenology - the relative amplitude of the order parameter between 
YBCO and LBCO (1381, the connection between charge order and magnetic instabilities 
(139 ■ 1411, the origin of finite correlations and the role of disorder in the chain layer (1421, 
the feedback on lattice dynamics (143 ■ 1451 - as a function of doping, temperature, and 
magnetic field. This tremendous amount of progress made it possible to shed new light 
on the role of charge order within the phase diagram and its connection to coexisting and 
neighbouring electronic orders and phases. The complete doping dependence of the salient 
properties of charge order in YBCO, and its comparison to La-based cuprates, were reported 
by Blanco-Canosa et al. (1321 and Mucker et al. { 1411 . 

One of the open questions is the connection between the charge order seen in high 
magnetic fields via NMR ( 129[|146 1 and the charge modulations seen in zero field via x-ray 
diffraction. A held-induced thermodynamic phase transition was detected in the sound 
velocity of underdoped YBCO (1471, suggesting that the high-held and low-held states are 
distinct, and in particular that the high-held charge order is two-dimensional in nature. 
NMR measurements, owing to their ability to probe both the high-held and the zero-held 
regimes, also established that charge order in these two regimes manifests itself as distinct 
phases, which furthermore coexist at low temperatures and high helds (1481. A recent x- 
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NdBCO: 

NdBa2Cu306+a; 


ray diffraction study in pulsed magnetic fields up to 28 T revealed that there is indeed a 
field-induced crossover from a short-ranged order to a long-ranged charge modulations with 
different period along the c axis (1311. 


3.2. Unifying real and reciprocal space: Bi2201 and Bi2212 

The discovery of CDWs in YBCO reinforced the idea that charge order might be a genuinely 
universal instability of the Cu02 planes, and the combination of real-space (STM) and 
reciprocal-space (neutron scattering, XRD, RXS) techniques had provided support for one 
and the same underlying general phenomenology in underdoped cuprates around 12 % hole 
doping. Despite such mounting evidence and the several unifying traits linking the various 
cuprate families, the real-space phenomenology of charge order in Bi-based compounds 
appeared very granular, at variance with the well-defined structures in reciprocal space as 
observed by scattering probes in La- and Y-based cuprates. The very different probing depth 

- few Angstroms for STM vs. hundreds of nanometer (and more) for scattering techniques 

- and a possible dichotomy between surface and bulk [as observed for charge-density waves 


along the Cu-O bond direction in LSCO (1491 and along the zone diagonal in Bi2201 


(150l], also contributed to a perceived disconnect between the domains (real space/surface 
vs. momentum space/bulk) and materials (Bi- vs. La- and Y-based cuprates) investigated 
by these two classes of experimental methods. 


Two recent works (49 1511 addressed this aspect by investigating charge order on the 


same materials using STM and RXS. Due to the complications in exposing atomically flat 
surfaces in any cuprate material other than Bi-based compounds, these two studies were 
performed on La-doped Bi2201 (in the doping range 0.11 <p< 0.14) ( [TST] ) and Bi2212 (in 
the doping range 0.07<p<0.13) (49l. In Bi2201, charge order was found with a wavevec- 


tor between 0.243 and 0.265 r.l.u. for decreasing doping, whereas in Bi2212 the ordering 
wavevector spans across a more extended range, from 0.25 to 0.31 r.l.u., a finding which was 
later corroborated by the detection of charge order also in optimally-doped Bi2212 (p~0.16, 
rc = 98K) ( |152[ ). Most importantly, both studies demonstrated that the LDOS modulations 
imaged by STM and the reflections measured by RXS in reciprocal space originate from the 
very same microscopic entity. This is shown in Fig. [7|r and b by the comparison between 
the RXS momentum scans (panel a) in underdoped Bi2201 (Tc = 15K, p~0.11), revealing 
a charge order peak emerging at low temperature at Q ~ (0.265,0,1/), and the Fourier- 
transformed differential conductance data (panel b) taken on the very same samples, which 
exhibit a peak at the same momentum value. A similar correspondence was found in Bi2212. 
The temperature dependence of the RXS charge order peak for three different Bi2201 dop¬ 
ing levels is reported in Fig.[^, showing a very gradual onset of the density modulations 
which occurs on a temperature scale which is proximate to the pseudogap temperature 
T* as determined from Knight shift measurements (1531. In Bi2212, RXS data also re¬ 
veal a rather slow temperature evolution of the peak intensity, which is accompanied by a 


drop below Tc- Although this drop is smaller than its analogue in YBCO (121 133 1361 


the competition with superconductivity is clearly demonstrated by temperature-dependent 
STM measurements that also reveal that the charge order is more pronounced for the un¬ 


occupied states (49l. Furthermore, in Bi2201 comparison between the RXS results and the 


Fermi surface measured using angle-resolved photoemission spectroscopy (ARPES) revealed 
a quantitative link between the observed charge order wavevector and the momentum vec¬ 
tor connecting the tips of the Fermi arcs, which in this case coincide with the so-called 
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Figure 7 

RXS and STM joint evidence of charge order in Bi2Sr2Cu06+5. (a) RXS scans along 
(iTOO) for an underdoped (Tc~15K) Bi2201 at low (blue) and high (red) temperature, showing 
the emergence of a broad charge order peak around i7~0.25 r.I.u.. (b) (HOO) linecut of the 
Fourier-transformed differential tunnelling conductance on freshly-cleaved surfaces of the same 
Bi2201 samples, indicating the equivalence of the charge modulations detected using STM. (c) 
Temperature dependence of the charge order peak intensity (from RXS) for three doping levels 
(Tc~15K, 22 K, and 30 K) in the under-to-optimal doping range, and correlation of the charge 
order onset with the pseudogap temperature T* (grey boxes), (d) Experimental (bottom) and 
calculated (top) Fermi surface for UD15K Bi2201, highlighting the Fermi arcs characterizing the 
pseudogap regime, and the wavevector (yellow connector) linking the arc-tips, or ‘hot-spots’, in 
agreement with the RXS experimental results. Readapted from Ref.|15l| 


‘hot-spots’ (see Fig. IT)- This correspondence suggests a possible link between the density 
modulations and the low-energy electronic structure, an element which is consistent with 
the doping evolution of the wavevector and which had been hinted by previous studies 
(46 1541. Recent theoretical works have discussed the possible special role played by the 
hot-spots in the context of a magnetically-driven charge order instability ( |13[|155}|158[ ), as 
well as the possibility that charge order might arise from a 2Q instability of the antinodal 
points within a pair-density-wave framework (159l. Additionally, we note that the RXS 
detection of charge order around (5~ 0.25 r.I.u. in Bi2201 confirms previously unpublished 
findings of a phonon anomaly at the same momentum location using inelastic x-ray scat¬ 
tering ( |160[ ) , which might be suggestive of a related anomaly in the electronic susceptibility 
or otherwise of a very strongly momentum-dependent electron-phonon coupling mechanism 
at work. 


www.annualreviews.org • Resonant x-ray scattering studies of charge order in cuprates 



















_ 200 
cu 

L_ 

13 

CO 

L_ 

<D 

CL 

0 100 


0 

0.35 


^ 0.30 
d 

lJ 

O 

o 0.25 
o 
CD 
> 

> 

CO 

^ 0.20 


0.151_ 1 1 I _ 

0 0.05 0.10 0.15 0.20 

Doping p (holes/Cu) 

Figure 8 

Charge order onset temperature and wavevector for all cuprate families, (a) Onset 
temperatures of charge order vs. hole doping, for all cuprate families. The shaded graphics in the 
background outline the boundaries of the antiferromagnetic (AF), superconducting (SC), and 
pseudogap (PG) regimes in YBCO. (b) Doping dependence of the charge order wavevector (the 
LNSCO and LSCO experimental points from Refs. |40||161| are calculated from the position of the 
spin ordering wavevector). Full symbols are from momentum-resolved probes (RXS, XRD, 
neutron scattering), while open symbols are from real-space methods (STM). Colored lines are 
guides-to-the-eye; the vertical dashed line marks the location of the doping p = 0.12. 
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3.3. Towards a universal phenomenology: charge order in Hgl201 and NCCO 
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In the wake of the results on Y- and Bi-based cuprates, evidence had been mounting that 
charge order could be a universal phenomenon in the cuprates. In 2012, Wu et al. detected 
charge order in the prototypical superconductor LSCO, around a doping of 12% (1491, a 
discovery which revealed that the low-temperature tetragonal structure in La-based cuprates 
is not essential for the appearance of stripe order (but perhaps relevant for stabilizing 
it). In this study, the comparison between resonant and non-resonant scattering data was 
interpreted in terms of a surface enhancement of the stripe order in LSCO. More recently, 
stripe order was found in the bulk of LSCO using hard x-ray diffraction (162 1631 and 
resonant x-ray scattering (164|. 

Having exhausted all other families, Hg-based cuprates remained the last hole-doped 
family to be investigated. In particular, HgBa 2 Cu 04 +i (Hgl201) represents the only se¬ 
ries of compounds with a pristine, tetragonal unit cell, therefore possessing the highest 
structural symmetry among all cuprates. The first indirect evidence of broken translational 
symmetry in Hgl201 came from the discovery of Fermi-surface reconstruction via high-field 
measurements of the Hall and Seebeck coefficients in underdoped Hgl201 (165 I, followed by 
the detection of quantum oscillations (1661. In 2014, Tabis et al. reported the very first evi¬ 
dence of charge order below 200 K in underdoped (p~ 0.09, = 72 K) Hgl201 from RXS and 

RIXS measurements ( |167| ). The charge order peak was found at Q ~ 0.28 r.l.u. - comparable 
to the values found in Bi2212 and intermediate between YBCO and La-based compounds 
- and required using resonant excitation, indicating that charge order is a rather subtle 
phenomenon in the Hgl201. This instance, combined with the fact that Hgl201 exhibits 
record-high Tc among single-layered cuprates, suggests a possible anticorrelation between 
charge order and superconductivity. Furthermore, the study revealed a direct connection 
between the onset of charge correlations and T**, “the temperature at which the Seebeck coef¬ 
ficient reaches its maximum value, and below which conventional Fermi-liquid planar charge 
transport is observed”. It also suggested a possible common origin of the main instabilities 
in the Cu 02 planes, namely “the possibility that the sequence of ordering tendencies (q — O 
order precedes charge order, which in turn precedes superconducting order) and the phase 
diagram as a whole are driven by short-range antiferromagnetic correlations”. Lastly, Tabis 
et al. successfully established the link between the charge order-induced reconstruction of 
the Fermi surface and the size of the electron pockets observed by quantum oscillations in 
YBCO ( 122[|123' I and Hgl201 (166), and were able to correctly predict the QO frequencies 
over an extended doping range. 

Figure provides a graphical compendium of all charge order observations in the 
hole-doped cuprates for what concerns the doping dependence of the onset temperature 
(Fig.|^) and of the wavevector (Fig.|^) as observed using both spatial- (open symbols) 
and momentum-resolved (full symbols) probes. The shaded phase diagram in Fig.j^ is rep¬ 
resentative of YBCO and includes the antiferromagnetic (AF) phase near zero doping and 
the superconducting (SC) dome at higher hole doping levels. All coloured lines are guides- 
to-the-eye, except for the Q = 2p line which interpolates the doping dependence of the 
ordering wavevector in LSCO below 12 % doping according to the picture of perfect stripes 
with 1/2 hole per unit cell along the charged rivers. Note that the points from Refs. 40{|161 


are from neutron scattering measurements of the antiferromagnetic peaks in LSCO, from 
which the charge order wavevector Qco was derived using the phenomenological relation 
Qco =2Qaf. It is clear from the phase diagram of Fig. how in all cases charge order 
emerges at the highest temperature scales around 12 % hole doping and spans a doping 


Hgl201: 

HgBa 2 Cu 04+5 

NCCO: 

Nd 2 - 2 ,Cea;Cu 04+5 
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range tentatively bound by two endpoint | |52[[5^|l69[|l^^ , thereby hinting at a possible 
intimate connection between this phenomenology and quantum criticality 1 171|172 1. Inter¬ 
estingly, the onset temperatures appear to scale inversely with the charge order amplitudes 
and correspondingly correlation lengths (which are maximum for La-based cuprates, inter¬ 
mediate for YBCO, and weak in the case of Bi-based compounds). At the same time, the 
extent of the suppression of Tc near 12 % doping is larger for the families exhibiting stronger 
charge order, consistent with the experimentally-observed competition between these two 
instabilities. Finally, we note that two main phenomenological differences exist between 
the stripy La-based cuprates and all other families: (i) in La-based compounds charge or¬ 
der and incommensurate antiferromagnetic spin order are simultaneously present, whereas 
they are mutually exclusive in the other compounds; (ii) the doping evolution of the charge 
order wavevector exhibits the opposite sign for La-cuprates with respect to the other fam¬ 
ilies, which are putatively compatible with a nesting scenario (see also guides-to-the-eye 
in Fig.[^). There are therefore strong indications for a common charge instability in all 
families of doped cuprates, but at the same time there are also stark differences between 
the manifestations of charge order in the conventional stripy compounds and in all other 
families. 

With charge order consistently detected across all flavours of hole-doped cuprates around 
12 % doping, the next frontier to be crossed was represented by the exploration of the 
electron-doped side of the phase diagram. Early insights came from inelastic x-ray scattering 


work (1731 reporting anomalies in the dispersion of optical phonons in Nd 2 -xCea:Cu 04+5 


(NCCO), followed by quantum oscillation studies providing evidence for a reconstructed 


small-pocket Fermi surface in the same compound (1741 and, more recently, by evidence 


from time-resolved studies of a broken-symmetry phase (175l. Recently, the first hints of 


charge order were reported by Lee et al. (1761 and Ishii et al. (1771, who used RIXS to 


map out the momentum-dependent structure of low-energy bosonic excitations in NCCO, 
for electron-doping values spanning the phase diagram from the antiferromagnetic to the 
superconducting phase. In both studies, RIXS measurements revealed the presence of spin 
waves associated to the antiferromagnetic order at low doping, which evolved into param- 
agnon excitations in the superconducting state (a; ~ 0.15). Around this doping, an addi¬ 
tional branch of a rapidly-dispersing, gapped excitation was found, which was interpreted as 
a particle-hole excitations in Ref. |177| and as a charge amplitude mode in Ref. |176| therefore 
suggesting the possible signature of a symmetry-broken quantum state over an extended 
temperature range beyond the superconducting phase. 

Direct evidence for charge order in NCCO was obtained by Da Silva Neto et al. ( |178| ) 
by means of RXS measurements in reciprocal space, which revealed a broad reflection at 
a wavevector Q ~ 0.25r.l.u., i.e. very proximate to the findings in hole-doped cuprates, 
once again pointing to a unified phenomenology and a possible universal instability of the 
Cu02 planes. The RXS scans for a superconducting NCCO sample (a; = 0.14) are shown 
in Fig.[^ and b as a function of photon energy and temperature, respectively. The charge 
order peaks appear to be rather broad, with a correlation length 25 — 35 A, and a very 
smooth evolution as a function of temperature and a finite peak amplitude seemingly per¬ 
sisting up to rather high temperatures (350-400 K), as reported in Fig.j^. Such a high onset 


^Note however that recently charge order has been observed in YBCO in the very underdoped 
regime, at p ~ 0.058 (Tc ~ 12.6 K), along the K reciprocal axis only, with wavevector Qco ~ 


0.337r.l.u. and onset temperature Tqo~65K (1681. 
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Figure 9 

Resonant scattering study of charge order in the electron-doped cuprate NCCO. (a, 6) 

Energy- (a) and temperature-dependent (6) RXS scans [along (ffOO)] of the charge order peak 
around if ~0.25 r.l.u. in superconducting Nd 2 - 2 :Cea;Cu 04_|_5 with 3; = 0.14. (c) Gradual 
temperature evolution of the charge order intensity, with a tentative onset temperature around 
350-400 K. Readapted from Ref. |178l 


temperature rules out a direct conuectiou of charge order with the pseudogap pheuomenol- 
ogy in electron-doped cuprates; instead, the charge order signal appears to partly correlate 
with the temperature evolution of the antiferromagnetic fluctuations in the same material 
( |179[ ), uncovering a possible connection between the charge and magnetic instabilities, as 
previously suggested for the hole-doped cuprates 1155 157 167 1801. Most importantly. 


and independently of the detailed temperature dependence, the detection of charge or¬ 
der in NCCO conclusively demonstrates that this phenomenon is truly universal in the 
cuprates, establishing a powerful and robust paradigm for the physics of the lightly-doped 
Cu 02 planes, which transcends the asymmetry inherent to the different orbital character of 
doped holes (largely on 0-2p states) vs. doped electrons (occupying the Cu-derived upper 
Hubbard band). 


3.4. What hides behind a peak: charge order patterns and symmetries 

As the case for the universality of charge order was receiving increasingly supporting evi¬ 
dence, the theoretical understanding of the origin of this phenomenon and of its interplay 
with coexisting instabilities, as well as the detailed exploration of the microscopic structure 
of the ordered state, all regained a central role in the context of the physics of copper-based 
high-temperature superconductors. |l3l [l55p5^pl|fl85t 

In a single-band model charge order would be a scalar field, which can be expressed as 
the periodic modulation of the occupation of a given electronic orbital as a function of the 
spatial coordinate. However, the low-energy electronic structure of the Cu02 involves three 
different orbitals: Cu-3d,j,2_y2 (at the Copper sites), 0-2p^ (at the horizontally-bridging 
Oxygen sites), and 0-2py (at the vertically-bridging Oxygen sites). Consequently, charge 
order can be expressed as a vector field with three components 1 156[|180 1: (i) a site-centered 
modulation, corresponding to an extra charge residing on the Cu-3d orbital (Fig. |10} :, top); 
(ii) an extended s’-wave bond-order, where the spatially-modulated density is on the 0-2p 
states, and the maxima along the x and y directions coincide (Fig.|10}:, middle); (hi) a d- 
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Figure 10 

Azimuthal angle-dependent charge scattering and symmetry of charge order in 

YBCO. (a,6) Azimuthal geometry for RXS experiments, which allows rotating the 
crystallographic axes of an angle a around the transferred momentum Q, which is consequently 
preserved in the frame of reference of the sample, (c) Different (but not orthogonal) symmetry 
representations for the charge ordered state in a 3-orbital system (Cu-3d, 0-2p^, and 0-2py), such 
as the cuprates, (d) The modulation of the RXS signal vs. azimuthal angle a is evident from the 
raw RXS scans across the charge order peak in underdoped YBa2Cu306.75, for both a and tt 
polarizations, (e) Azimuthal dependence of the ratio of the a and it charge order peak intensities, 
and comparison to the numerical prediction for different binary combinations of the symmetry 
terms introduced in (c), suggesting a prominent d-wave component. Readapted from Ref.|186| 


wave bond-order, where the charge modulation changes sign between x- and y-coordinated 
oxygen atoms, and the maxima are shifted by a half wavelength (Fig. |10[ ;, bottom). The 
notation, as well as the denomination of ‘symmetry terms’ for the components introduced 
above, follow from the resulting angular distribution of the phases within the CUO 4 plaque- 
tte; this is encoded via an internal momentum variable k (restricted to the first Brillouin 
zone) in addition to the lattice momentum Q, in the definition of the charge order param¬ 
eter: AcDw(k, Q) = ■ Ck-Q/z^ This definition enables a full decomposition of 

the orbital-dependent modulation of the electronic density as a linear combination of the 
symmetry components introduced above, which therefore serve as a basis set for the charge 
order parameter, so that the latter can be expressed as ( |156| ): 

AcDw(k, Q) = As -I- As> {cos kj:+cosky) + Ad{coska:—cosky) (15) 


where As, As', and A,; represent the magnitude of the s-, s’-, and d-wave terms. 

Theoretical predictions in the context of the t-J model 113 156 180| 188 1 and of the 
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Figure 11 

The symmetry of charge order in LBCO. (a) Azimuthal dependence of the RXS signal from 
a stripe-ordered LBCO sample, for the two incoming light polarizations a and tt, with theoretical 
fit profiles achieving best agreement in the case of predominant s' symmetry. (6) Ratio I-jr jla of 
the intensities for the two polarizations, and corresponding best fit profile, (c) Diagram of the 
relative weight of s'- and d-wave components as a function of the sign and magnitude of the 
in-plane (iy) to the out-of-plane {t±) x-ray transition amplitudes; the red marker and bar indicates 
the parameter range yielding the best agreement with the data. Readapted from Ref.|187| 


spin-fermion model (155 1581 concluded that a d-wave pattern of electronic charges is 
energetically favoured over the other terms. In order to test these theoretical scenarios, an 
alternative RXS scheme has been recently devised 1186 187[ ) and applied to the study of 
the local symmetry of charge order in Bi2201 and YBCO |l86[ ), and YBCO and LBCO 
(1871. This approach relies on the definition of a local form factor fpq (which is a tensorial 
quantity, see Eq.[^ capable of discriminating between the different symmetry components 
of charge order. Once the magnitudes of the different components are built into the local 
form factor [fpq —>■ fpq (As, As', Ad)], the scattering yield can be written as: 


(As, As', Ad) 


(16) 


where and represent the polarization vectors for incoming and outgoing photons, 
respectively, while Q* is the ordering wavevector. At this point, a single measurement of 
the RXS intensity will not suffice to resolve the As, As', Ad terms. This issue is overcome by 
performing the RXS measurement as a function of the sample rotation around the azimuthal 
axis (see schematic of experimental geometry in Fig. 10 r and b), which is collinear with the 
ordering wavevector. This procedure allows modulating the projection of the form factor 
onto the light polarization as a function of the azimuthal angle a. Following this approach, 
the RXS scans can be measured for a range of azimuthal values (see Fig. |10{ l for the case 
of YBCO, from Ref. |186p . The resulting RXS intensities, shown in Fig. |10[ ; for the case of 
YBCO at the Cu-Ls edge (1861 and in Fig. and b for the case of LBCO at the 0-K 
edge (1871, can be htted to the theoretical model of Eq|16| to evaluate the magnitudes 
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of the symmetry components (which are treated as variational parameters). From these 
studies, it was concluded that YBCO possesses a prominent d-wave symmetry (Fig.|10^ and 
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Ref. 1861, while the stripe order in LBCO is best described by a s’-wave pattern (Fig. 
and Ref. 187). The determination of the charge order symmetry in Bi2201 was instead not 
conclusive using RXS (186l. While these early RXS results are poised to stimulate further 
work to identify and classify these types of symmetries in momentum space, in a recent 
STM study of Bi2212 and Na-CCOC, Fujita et al. (|189|) have successfully implemented the 


decomposition of Eq.[^ by analyzing the reciprocal space representation of the different 
local symmetries in real space. In particular, these authors first noted how a d-wave [s’- 
wave] symmetry suppresses the Fourier amplitudes of the (±<5, 0) and (0, ±Q) [(±1 ± Q, 0), 
(0, ±1 ± Q), (±1,±(3), and (±(3,±1)] peaks, and subsequently were able to resolve the 
extinction of the inner charge order peaks in reciprocal space starting from the spatially- 
resolved, partial orbital occupation of the Ox, and Oy sites. The detailed analysis of the 
STM conductance maps finally revealed that charge order in Bi2212 and Na-CCOC also 
possess a predominant d-wave form factor. 

Another key aspect regarding the microscopic description of charge order, and one that 
has been debated for long time ( 7| |12[[T83l|190 1941, is whether charge order has a checker¬ 
board (bidirectional) or stripe (unidirectional) character. This problem found an early 
answer in the La-based cuprates, thanks to the coexistence of spin and charge order with 
a precise wavevector relation that rules out a checkerboard state. However, magnetic and 
charge order tend to avoid each other in the phase diagram of all other cuprates, a cir¬ 
cumstance which, combined with the typical observation of orthogonal [(Q, 0) and (0, Q)] 
charge order reflections, has hindered a conclusive resolution of the checkerboard/stripe 
dichotomy. Furthermore, evidence in real space using local probes (STM) has been tradi¬ 
tionally hampered by the disorder characterizing Bi-based compounds, known to blur the 
distinction between native uni- and bi-directional ordering instabilities 1 192[|193 1; however, 
recent advancements in the analysis of STM datasets have enabled the visualization of a 
predominant unidirectional character of electronic modulations in Bi2212 ( |195| ). 

The first indications of unidirectional charge order in underdoped YBCO came from 
Blackburn et al. (139) and Blanco-Canosa et al. (1401, showing experimental evidence for 
very unequal amplitudes between the charge order peak along the b axis (strong) peak and 
along the a axis (weak) in YBCO Ortho-II, a direct signature of a tendency to unidirec¬ 
tional behavior in the charge order parameter. A recent attempt to assess the character 
of density modulations was put forth by resolving the full reciprocal-space structure of the 
charge order peaks in YBCO using RXS (1961. The same geometry was used as outlined in 
Fig. 10 1 and b, which effectively allows to slice through the ordering peak in the iQx,Qy) 
plane, along different directions in reciprocal space. In this study, the linewidth from the 
RXS scans (see again Fig. 101) was extracted in order to reconstruct the two-dimensional 
peak shape as shown in Fig.|12^ for a YBa2Cu306,5i sample. The elongation of the charge 
order peaks located along the (100) direction (Qa) and the (010) direction (Q;,) suggests a 
locking between the direction of long correlation (narrow peak linewidth) and the wavevec¬ 
tor, occurring for the YBa 2 Cu 306 . 5 i and YBa 2 Cu 306.67 samples - whereas the charge order 
peaks in YBa2Cu306.75 exhibits the same elongation. This locking mechanism suggests a 
breaking of four-fold rotational {D 4 h) symmetry which is independent for the two ordering 
components, and is therefore incompatible with a checkerboard state - the latter being an 
equal superposition of density modulations along (100) and (010), imposes the same peak 
structure for the charge order peaks Qa and Q;,. 
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Figure 12 

Stripe vs. checkerboard symmetry of charge modulations in YBCO. (a) 

Two-dimensional shape of the charge order peaks at (0.31, 0,1.5) (red ellipse) and (0, 0.31,1.5) 
(blue ellipse) for a YBa 2 Cu 306 . 5 i sample, as determined from fitting the azimuthal-dependent 
peak widths (red and blue markers). Bottom-right inset: schematic representation of the original 
peak shapes (full ellipses) and their 90° rotated versions (hollow ellipses). (6) 
Temperature-dependent longitudinal (green) and transverse (yellow) correlation lengths at 
(0.31, 0,1.5) (left) and (0, 0.31,1.5) (right), showing a clear anisotropy in the evolution across Tc. 
(c) Illustration of the anisotropy reported in (6), showing how the onset of the superconducting 
phase reduces the density-density correlations preferentially across the stripes. Readapted from 

Ref.im 


The momentum structure of the charge order peaks already rules out a checkerboard 
state in favour of a unidirectional (stripy) instability (where stripes can be segregated or 
overlapping and still generate the same structure in reciprocal space). Further support is 
provided by the temperature dependence of the longitudinal (^y) and transverse (^x) corre¬ 
lation lengths, which can be derived as the inverse of the peak linewidth along the direction 
parallel and perpendicular to the wavevector, respectively, and are shown in Fig. |12[ 3. The 
temperature evolution of the correlation lengths suggests a larger drop, below the supercon¬ 
ducting transition temperature Tc, for the longitudinal (i.e., across the stripes) correlations 
(see cartoon in Fig. 12:), again leading to a preferential, q-dependent locking of the density- 


density correlations and to a breaking of fourfold rotational symmetry. 

A similar tendency to a unidirectional character of the charge modulations has also been 
identified in the different orbital symmetry of the Qa and Qi, peaks, measured with RXS 
(|187[). This in-plane anisotropy of the short-range charge modulations in YBCO, also de¬ 


tected in NMR measurements (1481, have been argued to cause the large in-plane anisotropy 
of the Nernst signal seen in YBCO near p = 0.12 (197l, since the Nernst anisotropy grows 
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upon cooling in tandem with the growth in modulation amplitude (1981. 


4. FUTURE PROSPECTS AND NEW CHALLENGES 

Despite the recent flurry of experimental findings and theoretical insights on charge or¬ 
der in the cuprates, there is still a lot to learn and to understand. First and foremost, 
the mechanism driving the doped holes into breaking translational symmetry has not been 
conclusively pinned down. In particular, it is unclear whether the cuprates exhibit Peierls 


physics similar to other low-dimensional compounds (1991 or whether a new and uncon¬ 


ventional mechanism is at play. Several studies have been performed since the discovery of 
stripe order, aimed at elucidating the electron-lattice interplay and its relevance for charge 


order instabilities (200 ■ 2051. Recently, the work by Le Tacon et al. (1441 revealed a strong 


and sharp (in momentum) softening of the low-energy acoustic and optical phonons at the 


charge order wavevector in underdoped YBCO [see Fig. 13 1 - similar observations were also 


made for NCCO (1731 and Bi2201 (160l]. This result per se already provides direct evi¬ 


dence of a pronounced electron-lattice coupling mechanism, while the partial yet incomplete 
softening (with the frequency remaining nonzero) further clarifies that the charge order in 
cuprates is not due to a phonon mode freezing into a static lattice distortion, consistent 
with the short-ranged nature of charge correlations at zero magnetic field. The charge order 
peak remains confined to the quasi-elastic line, which bears the typical signatures of a ‘cen¬ 
tral peak’, which is characteristic of the presence of ordered nanodomains which gradually 
coalesce into a state with longer correlations. However, the most puzzling finding emerges 
from the temperature dependence of the softening effect, which is found to be enhanced in 
the superconducting phase, where charge order is paradoxically weakened. This behavior 
reflects a strong superconductivity-induced phonon anomaly, occurring at the wavevector of 
the charge instability, thus exposing a very complex intertwining between these competing 
phases and might deserve further explorations for a recent theoretical discussion of this 


phenomenon see Ref. 2061. 


Time-resolved experimental methods, which probe the recovery dynamics of the sys¬ 
tem following an ultrafast perturbation, have rapidly emerged as an alternative tool to 
gain further insights on correlated materials. Depending on the nature of the excitation, 
pump-probe techniques can probe both the microscopic coupling between different degrees 
of freedom near equilibrium (i.e., in the linear response regime attained at low excitation 
intensities) ( 208f 2111 or explore completely new physics in the strongly perturbative regime 
(at high excitation intensities) where photoinduced phases can be accessed that are often in¬ 


accessible at equilibrium (2121. The field of time-resolved scattering and diffraction has seen 


rapid advancements since its inception at the end of the 1990’s (213l. New avenues in the 


x-ray study of ordering phenomena out-of-equilibrium have been enabled by the recent de¬ 
velopment of free-electron-laser (FEL) sources, which generate bright ultrashort light pulses 
over a broad range of energies from the extreme ultraviolet to the hard x-rays, thereby cre¬ 
ating new opportunities for the investigation of structural and electron dynamics with high 
time resolution. The spectral control of the photoexcitation process is an additional experi¬ 
mental parameter, enabling the selective perturbation of specific degrees of freedom via the 
nonlinear coupling between light and lattice ( 214||215' l. When this framework is combined 
with photon energy tunability, resonant scattering can be used to probe the recovery of 
electronic orders that are brought out-of-equilibrium by an ultrafast pump (216 2171. In 


the context of charge order in the cuprates, these measurements were performed for the 
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Figure 13 

Giant phonon anomaly and charge order melting from high-resolution frequency- and 
time-domain x-ray spectroscopy, (a) Momentum and temperature dependence of the 
frequency of low-energy tranverse acoustic and optical phonons in underdoped YBCO (Tc = 61K), 
showing a sharp softening occurring at the charge order wavevector Qco = (0) 0.31, 6.5) in the 
superconducting phase (T < Tc). Readapted from Ref. |144| (ft) Charge order peak in underdoped 
YBCO (rc = 61K), acquired using ultrafast soft x-ray pulses from the Linear Coherent Light 
Source, (c) Corresponding time-resolved dynamics of the charge order peak amplitude following 
impulsive photoexcitation (t = 0) of the apical oxygen mode (A~15/rm) using 400 fs pulses. 
Readapted from Ref.|207| 


first time by Forst et al. (2071, using ultrafast THz light pulses to photoexcite a particular 
lattice vibration (the apical oxygen mode) in underdoped YBCO and soft x-ray FEL pulses 
(at the Cu-La) to track the evolution of the charge order peak (Fig.[T^) as a function of 
time. The photoexcitation process strengthens the superconducting state by promoting co¬ 
herent interlayer transport and the charge order amplitude is correspondingly weakened by 
a factor two (see Fig. 13:), thus providing new insights on the competition between charge 
order and superconductivity in a regime where phase coherence is artificially enhanced with 
light ( 218||219] ) [similar results were also reported in LBCO (220 (]. More in general, the new 
FEL capabilities are poised to considerably extend the domain of application of resonant 
x-ray methods toward the study of nanoscale ordering phenomena at unprecedented length- 
and time-scales. 

Another new frontier for the use of high-brightness x-ray sources is in the context of 
high-magnetic field studies. The use of high fields within different experimental techniques, 
such as quantum oscillations or NMR, has been transformative for our understanding of 
high-temperature superconductors, enabling unprecedented insights into the nature and in¬ 
terplay of competing orders in these complex materials. It was then natural to envisage the 
extension of high-held capabilities to scattering and diffraction experiments and, notwith¬ 
standing the complications caused by the requirement for optical access for incoming and 
outgoing photon beams, several efforts have been successfully carried forward in this direc¬ 
tion in recent years. Chang et al. ( |133[ ) used a superconducting cryomagnet specihcally 
designed for low-angle forward-scattering measurements in transmission geometry and us¬ 
ing very penetrating high energy photons (lOOkeV), which enabled the hrst observation of 
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magnetic field-induced enhancement of the charge order signal. 

Very recently, an innovative experimental scheme was introduced by Gerber et al. 
that is based on the use of a high-field (28 T) split-pair pulsed magnet synchronized with 
the ultra-bright and ultrafast x-ray pulses (with photon energy of 8.8 keV) of the Linear 
Coherent Light Source. This approach, leveraging the high single-pulse photon flux of the 
FEL beam, enabled the acquisition of two-dimensional single-shot diffraction patterns with 
sufficient photon counts in spite of the very low sampling frequency (with one spectrum 
every 2 to 25 minutes, due to the recovery time of the pulsed magnetic field apparatus). 
The momentum-space maps of charge order in YBCO, from 0 to 25 T, reveal an ostensible 
evolution in the momentum structure of the charge order peak as a function of both the 
in-plane [K) and out-of-plane (L) wavevectors. In particular, the out-of-plane character of 
the charge order peak changes from a very broad, elongated structure at zero field to a well- 
defined peak (with ~170A correlation length) located around L = l. This is accompanied 
by a sharpening of the in-plane peak linewidth, with a corresponding increase of charge cor¬ 
relations from ~60 to ~180 A, and a concomitant enhancement of the diffracted intensity. 
These results disclose the complexity of charge order and its rich and unconventional phe¬ 
nomenology as a function of doping as well as magnetic field. In particular, this study not 
only uncovers the momentum-resolved crossover between the low- and high-field regimes, 
previously investigated using NMR (1481, but also clarifies that a full reconstruction of the 
Fermi surface as seen by quantum oscillations only happens in high field when long range 
order sets in, thereby elucidating the absence of folding and small pockets in the ARPES 
measurements (which can only be performed in zero field) on YBCO ( 221[|222 1. 

To conclude, we would like to emphasize that the further development of these novel 
approaches and methodologies will enable the exploration of completely new dimensions. 
This will bring us many more surprises and much deeper insights in the study of symmetry 
breaking instabilities in cuprates, and will lead to a fuller understanding of these phenomena 
and of their intimate interplay with high-temperature superconductivity. 


SUMMARY POINTS 

1. Technical advancements: The challenge posed by the detection and characteriza¬ 
tion of charge-density waves in the context of the rich phenomenology of the doped 
Cu 02 planes has propelled tremendous advancements in the field of soft x-ray scat¬ 
tering methods. Presently, RXS beamlines are operational or under construction 
at several facilities worldwide: ALS, APS, BESSY, CLS, DESY, Diamond, ESRE, 
NSLS-II, NSRRC, SLS, SOLEIL, Spring-8, SSRL. 

2. Universality: To-date, evidence of charge order has been found for all hole-doped 
cuprate families, as well as in Nd-based electron-doped compounds, using a variety 
of complementary experimental probes. 

3. Resolution: Resonant x-ray scattering has been successfully used to detect charge 
density modulations with a spatial coherence as short as 15-20 A, thus extending 
to momentum space capabilities that were previously accessible only to spatially- 
resolved techniques (STM). 

4. Symmetry: The richness of the RXS information, and the multiple control pa¬ 
rameters available - polarization, photon energy, sample orientation - have opened 
up the possibility to explore new microscopic aspects of the ordered state - such as 
its dimensionality (ID vs. 2D) and its local symmetry (s-, s’-, and d-wave). 
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FUTURE DIRECTIONS 


1. Unveiling the driving force behind charge order in cuprates. 

2. Find new approaches to modulate and, in general, control charge order. 

3. Elucidate the relative role of Mott physics and the proximity to quantum-critical 
behavior for charge order, as well as superconductivity, in cuprates. 

4. Probe the evolution of charge order and its interplay with superconductivity under 
extreme conditions - high magnetic fields, high pressures, ultrafast optical pumping. 
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